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Abstract 

The basin of attraction of a stable equilibrium point is investigated for 
a dynamical system (W97) that has been used to model transition to 
turbulence in shear flows. The basin boundary contains a linearly unstable 
equilibrium point Xib which, in the self-sustaining scenario, plays a role 
in mediating the transition in that transition orbits cluster around its 
unstable manifold. We find for W97, however, that this role is played not 
by Xib but rather by a periodic orbit also lying on the basin boundary. 
Moreover, it appears via numerical computations that all orbits beginning 
near Xw relaminarize. We offer evidence that this is due to the exquisite 
narrowness of the complementary region to the basin of attraction in 
the part of phase space near Xw- This further leads to a proposal for 
understanding the 'edge of chaos' in terms of more familiar invariant sets 
of a dynamical system. 
MSG numbers:76D05, 76F20 



1 Introduction 

Experimentally, laminar shear flows undergo transition to turbulence when the 
relevant parameter, the Reynolds number R, exceeds a critical value Rc. Math- 
ematically, when the Navier-Stokes equations are linearized about the laminar 
flow, the expected passage from stability to instability at Rc is not found. This 
is the familiar conundrum that linear theory fails to predict the critical value 
Rc (cf., for example, the introductory remarks in 5 for a fuller discussion). A 
resolution of this conundrum is that the stable, laminar point O possesses a 
basin of attraction B whose boundary dB passes increasingly close to O with 
increasing R, so that perturbations that may be small by laboratory standards 
are large enough to transgress dB for sufficiently large values of R. 

This idea of describing the problem of shear flows in the language of dynamical- 
systems theory is an attractive one but has limitations when the model systems 
are confined to very low dimensions. For example, the debate whether turbu- 
lence is transient or not (|4]) can hardly be joined in this context, where dB 
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presents a clear boundary between transient and permanent departures from 
the laminar flow. The present work is nevertheless confined to models of very 
low dimensions. The motivations are (1) the impression that the nature of the 
basin of attraction is an important element in the theory, (2) the observation 
that very little is known about it and (3) the conviction that it would be a good 
idea to understand the basin and its boundary in low-dimensional systems be- 
fore proceeding to high-dimensional systems. In this paper we have considered 
the four-dimension model W97 (as described in [9J) and minor modifications of 
it. 

Some results that may be relevant to higher-dimensional models and to the 
Navier-Stokes (NS) equations, discussed in more detail in ij5l include the struc- 
ture of the boundary, which implies that the functional form of a perturbation 
may be as important as its size; that the vicinity of X/f, may not be a good 
place to seek transition; and that the tendency of the complementary region to 
the basin of attraction to extreme narrowness in some parts of phase space may 
help to explain the 'edge of chaos' (fZ]) in terms of the more familiar invariant 
sets of dynamical-systems theory. 

The plan of the paper is as follows. We describe the mathematical setting 
in 3H In iJ3]we present WalefFe's model together with diagrams of the boundary 
of the basin of attraction indicating the periodic orbit that lies on that bound- 
ary, and the relaminarization of orbits starting near X/^. In |]3]we indicate a 
resolution of this relaminarization in terms of the folded structure of the basin 
boundary. The concluding section, [JSJ is devoted to drawing from this resolu- 
tion a proposed interpretation of the edge of chaos and to brief remarks on the 
results of this paper and their implications for further study. 

2 Mathematical Setting 

The Navier-Stokes (NS) equations possess a number of very simple solutions rep- 
resenting laminar shear flows (plane Couette and Poiseuille flow, pipe flow, etc.). 
When these partial-differential equations are modeled by a finite-dimensional 
system, the laminar flow can modeled by an equilibrium point of that system, 
which we'll take to be the origin of coordinates O. Almost all such finite- 
dimensional systems that have been studied take the form 

i = Aa; + 6(x),xei?" (1) 
satisfying certain conditions: 

1. A is a non- normal, stable matrix. 

2. h{x) is quadratic in x and X)j=i ^j^ji^) = 0- 

This structure can be inferred by a Galerkin projection of the NS equations onto 
n basis vectors, while taking mild liberties with the boundary conditions. 

The stability condition on A implies that its eigenvalues lie in the left half- 
plane and therefore that the origin O is asymptotically stable. The basin of 
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attraction B of an asymptotically stable equilibrium point is the set with the 
property that any orbit through a point of B tends to the equilibrium point 
as t +00. It is an open set invariant under the flow, in the sense that 
any solution beginning in B remains in B on its maximal interval of existence 
(a, 00 j3- It's boundary dB (if it has one) is likewise invariant in the same sense. 
The latter is typically of relative measure zero so orbits that lie on dB are 
rare but important since they lie just at the transition from the laminar flow 
toward something "more interesting." In particular, the threshold amplitude 
for transition, which has been a subject of some interest ([10], 12, [J), is the 
minimum distance from the origin to the basin boundary. We denote such a 
threshold point by T. 



Figure 1: This cartoon shows Waleffe's picture; a part of the boundary dB of the 
basin of attraction of the point O is shown. This part of the boundary coincides with 
the stable manifold of the equilibrium point U. Orbits starting near O but on dB are 
attracted toward U. Those starting near O but slightly above dB are attracted first 
toward U but are then captured by the unstable manifold of the latter and carried 
away (toward turbulence). An example is the orbit through the threshold point T. 



3 The model W97 

This is a four-dimensional model which may be written 



^In some examples a 



Xi = 


—SriXi 


+ CTlxl - 


' Cr2X2X3, 


(2) 


X2 = 


-5r2X2 


+ cr2.T3 4 


- a2XiX3 - CT4X4, 


(3) 


±3 = 


-Sr^xs 


+ (^3x1, 




(4) 


±4 = 


-(o-i + 


Sr4) X4 -\ 


~ X4 {(74X2 - <7lXl - CT3X3) . 


(5) 


les a = 


—00, though this is 


not inevitable. 
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Here 6 ~ 1/R where R is the Reynolds number and the eight constants ri 
through (74 are all positive H. There are standard values for these constants (cf. 
[5]) that we use in this section. 

This system conforms to the rules ([T]) and ^ of model-building. It pos- 
sesses the symmetry 5* = diag(l, 1, 1, —1) so a solution x{t) has a companion 
solution x{t) obtained by reversing the sign of X4{t) so the plane 0:4 = is an 
invariant plane. For these reasons there is no loss of generality in considering 
only solutions for which 0:4 (t) > 0. It is not difficult to show that the invariant 
plane — lies entirely in B, the basin of attraction of the origin. 

When one seeks equilibrium solutions other than the origin, they are found 
in pairs provided 5 < 6sn or R > i?srill. The one lying closer in norm to the 
origin is called the "lower branch" equilibrium solution Xib, the one lying far- 
ther away is called the "upper branch" Xub (see Figure [2]) . The lower-branch 
solution is unstable, with a one-dimensional unstable manifold and a three di- 
mensional stable manifold; the upper-branch solution may be stable or unstable, 
depending on the choices of R and of the other parameters. For the values of 
R considered in this paper and for the standard values of the other parame- 
ters, Xub is asympotically stable. These equilibrium solutions are illustrated in 
Figure 




Upper branch 



Lower branch 



Figure 2: This is the diagram for equilibrium figures in Waleffe's model. The lower 
edge, where the norm ||X|| — 0, represents the laminar solution. The "lower branch" 
of further equilibrium solutions does not approach the laminar solution as ii ^ 00 
but instead tends to the indicated asymptote. However, as shown in [2], the curve of 
threshold values does approach the laminar solution, like . 



Figure ((T]) suggests that orbits starting on dB tend toward Xib, which me- 
diates the transition. This will be so if the the stable manifold of Xib coincides 
with dB. This seems plausible (and has been found to be the case for some 
other models) but is by no means inevitable: the basin boundary and the stable 



Hn an earlier model, W95, tri =0 ([§]). 

•^The value of Rg„ of course depends on the other parameters. With the standard choices 
for W97, Rsn = 104.84. 
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manifold of Xn, are both invariant sets for the system ^ but they are defined 
in different manners and need not be identical. In fact, we find that they are 
not identical for the model W97. The stable manifold of Xn, is a proper subset 
of dB, so a point on dB lying far enough from Xn, has a different evolution. 

For W97 the threshold point T has been located by Cossu ([2]). We find, by 
following orbits starting near T, that they are attracted not toward Xit but to 
a periodic orbit P also lying on dB. We have thus far carried out calculations 
for R = 145 and R = 190 and we display only those for R = 190 (those for 
R ~ 145 are similar). We exploit Cossu's calculations to find refinements of 
the threshold values T: by repeated bisection we obtain a pair of values, Xo 
and Xi lying respectively just outside and just inside the basin of attraction B 
and within a short distance e of one another. It follows that there is a point 
T on dB within e of either of them. We then find the structure of the basin 
boundary by calculating slices of it by various hyperplanes (Figures [3] and [5|) . 
We also obtain the orbit through T, finding the following. On taking Xo as initial 
data, we find that after a transient of a few hundred units of time, the orbit 
is essentially periodic for many thousands of units of time, eventually spiraling 
into the stable, outside point X^b', if Xi is taken as the initial point a similar 
evolution is found except that at the end, the orbit tends to the laminar point 
O. Neither of these orbits comes close to Xn,. These remarks are illustrated in 
Figures ©, © and ©. 




Figure 3: A slice through the basin boundary made by the hyperplane X2 = T2, X3 = 
T3, where Ti,i — 1, ... ,4 are the components of the threshold point T. The parts of 
this region lying in B are so marked; the remainder lies outside. An orbit starting very 
near T spends a very long time hovering near a periodic orbit, indicated in projection 
onto this hyperplane (the transient leading from T to this orbit is not shown). The 
complementary region to the basin becomes very narrow near T but is easily resolvable 
numerically for these values of the parameters. . 

We illustrate the results in the form of slices formed when certain two- 
dimensional hyperplanes intersect dB. Also seen in these diagrams are projec- 
tions of the periodic orbit - and other features - onto the hyperplanes in ques- 
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Figure 4: A global view of the preceding figure, showing more of the nature of the 
basin boundary, in the same slice as in Figure Q. 




Figure 5: Like Figure Q, but a different slice. 

tion. The projections are indicated by placing single quotation marks around 
their labels. 

While the self-sustaining process may be only slightly modified by replacing 
the equilibrium point Xn, with the periodic orbit P as the mediator of transition, 
the question of what role Xib plays in the dynamics now arises. We next turn 
to this. 

4 Orbits starting near the point 

The three-dimensional stable manifold of Xu, coincides locally with dB. Denote 
by ft a unit vector transverse to dB (for example, ft could have the unstable 
direction at Xu,). Then if for a scalar v we choose initial data 



x(0) = Xib + v^t 



(6) 



for a small value of we expect the orbit to depart from Xu, along its unstable 
manifold. We anticipate that for one sign of v the orbit will lie inside B and 
for the other outside, and therefore that for one sign the orbit will tend to the 
origin and for the other will remain permanently outside the basin boundary, 
presumably tending for large t to the stable equilibrium point at Xub- 

Instead, we find that all orbits tend to O, apparently echoing the persistent 
relaminarization found in other models (^). This is illustrated in Figure ([6]). 




200 400 600 800 1000 1200 200 400 600 800 1000 1200 

t t 



Figure 6: The model is W97 with the canonical values for the constants ri through 0-4 
and R — 190. The norms of orbits are shown on the interval [0,1200] and v — ±10~^^. 
They differ in the nature of the orbits and in the time for relaminarization to occur. 

This violation of expectations could be explained by the following conjecture 
for the system W97 near the unstable equilibrium point Xu,: X/f, indeed lies on 
dB but, near the part of dB on which it lies, there is a second leaf of dB 
exquisitely close to the first, and it's only for initial data in the narrow gap 
between the two leaves that orbits remain bounded away from the origin Q If 
this "gap conjecture" is correct for W97 with the standard choice of parameters, 
the space between the two leaves is too small to be detected in the usual double 
precision arithmetic. The strategy we adopt below for testing this conjecture is 
the following. 

If we put all the positive constants in W97 equal to unity we find qualitative 
similarity to the case with standard values. In particular, all perturbations of 
Xib relaminarize. The strategy will be to take all coefficients equal to unity 
with the exception of cti. An asymptotic analysis like that of [5] shows that, for 
small 5=1 / R., the lower and upper branch equilibrium points are 

Xif, (- 0-1(5, gi 6) 

1 -p (7-^ ± ~r (J ^ 

and 

Xub^{^l + 26,S'/',S'/',S'^'). 



*That narrow gaps are plausible for these systems may be seen by examining Figure ((Sjl 
near the point T. In that case the gap, while narrow, is still easily detectable numerically. 
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For small values of cti, we find numerically that there is a large and easily 
detectable gap near Xn,. We then consider successively larger values of cti to see 
if this gap gets successively narrower and ultimately becomes undetectable. In 
Figures [7] and [51 we show only slices of the basin boundary with the hyperplane 
X2 = Xib2 , X3 = Xib^ since these seem to reveal the gap most clearly. The value 
of R is held fixed at 15. 




Figure 7: The ordinate is X4,. For each value of ai shown, there is a progressively 
narrower gap to the right of Xn,. Orbits starting above the narrower gap require longer 
to relaminarize than those starting below. It looks on this scale as if the gap closes 
when (Ji = 0.185, with Xit falling on a critical point. A close up of the region near 
Xib would show that this has not yet occurred, but see Figure [S] 



These diagrams seem to confirm the gap conjecture but reveal a further, 
unexpected feature: there appears to be a topological change in the nature of 
the basin boundary for precisely the parameter value at which the gap becomes 
suddenly undetectable. We comment further on this in the discussion section 
below. 



5 Discussion 



In the self-sustaining scenario for the NS equations, the streaky flow is repre- 
sented by a steady-state solution (or by a traveling-wave solution: steady in 
a moving frame). These have been sought and studied in some detail in the 
context of the NS equations ([B],[TT]). The analog in W97 is the lower-branch 
equilibrium point Xn,. For standard values of the parameters of W97 we find 
that all orbits beginning near Xn, relaminarize, whereas in a different part of 
phase space - nearer to the threshold point - there is a large region wherein or- 
bits are permanently bounded away from the origin: the low-dimensional analog 
of persistent turbulence. Looking for this large region by beginning near Xn, 
would be counter-productive, and it is possible that the analogous conclusion 
holds for numerical studies of the NS equations (cf. the recent exploration of 
'relative periodic orbits' [3 in the NS context). 
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Figure 8: The large dots represent Xw- These are closeup views of the part of phase 
space near Xit- The region marked C is the complement of B. The first and third 
of these diagrams are calculated but the central diagram is conjectural. The curve 
marked 'edge' in the third diagram divides long-time relaminarizations from short- 
time relaminarizations: no points of the complementary region are detectable in the 
region of phase space immediately adjacent to this edge. 



The apparent collapse of the complementary set to the basin of attraction, 
as depicted in Figure is reminiscent of the 'edge of chaos' as described in 
(H) in that the line marked 'edge' is determined in a similar manner: there is 
a sharp difference in relaminarization time for points just above and just below 
this line. It suggests a picture like that of Figure wherein the 'edged in 
fact consists of two leaves of the basin boundary so close together as to be 
numerically indistinguishable. 

The diagram shown envisions a region of phase space close to Xib but that 
is because we have concentrated on this point. There are undoubtedly parts 
of phase space farther from Xn, where such narrow gaps occur, with behavior 
that is qualitatively similar but for which the transient relaminarizations may 
be quantitatively quite different. 

We tentatively propose the structure shown in this diagram as a building 
block for the edge of chaos as described elsewhere for other models (cf. [7]). 
This interpretation would require a repeated folding and refolding of the basin 
boundary of exquisite tightness, each folding modeled by such a building block. 
That intricate folds are possible for the basin boundaries of this and related 
systems is plausible not only on the basis of the pictures shown here but also 
experience with other model systems (e.g., [?])• It has been argued that this 
'edge' is itself a new kind of invariant set but it seems difficult to make this 
statement mathematically precise. In our proposed picture, the edge is from any 
numerical standpoint effectively invariant since orbits are as close as numerically 
possible to the basin boundary, which is indeed invariant. 

^There is no evidence of chaos in our calculations so we refer to it simply as the edge. 
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to Obstacle 



to ^ub Complement of B 




Figure 9: A conjecture regarding tiie 'edge.' The distance e separating the top and 
bottom leaves of the basin boundary is so small as to escape numerical detection. If 
a point could be located between these leaves, it would ultimately be captured by the 
stable equilibrium point Xub and therefore be permanently bounded away from the 
origin. Attempts to locate such a point are frustrated by the narrowness of the gap, 
and will result instead in initial points indicated by the asterisks above and below. 
The orbit starting above is complicated but, since it is in B, ultimately decays to O. 
The orbit starting below decays more directly, and more quickly, to O. 

A number of issues for further study come to mind in view of the outcomes 
of the work reported here. We hst a few: 

• What is the nature of the singularity indicated in Figure (O? It is possible 
in principle that the singularity is an artifact of viewing the surface dB 
through hyperplane slices and therefore does not represent any singularity 
of dB. Against this is the (apparent) fact that the equilibrium point Xn, 
of the vector field f{x) — Ax + b{x) lies at the singularity - an unlikely 
coincidence. 

• Can we exploit the picture presented here to locate regions of phase space 
where there is no relaminarization? One could do that in the context of the 
present model by following a long-time, relaminarization orbit, choosing a 
point on it which is safely far from both Xn, and from O, and conducting 
a search for dB near this point. 

• The structure of the basin boundaries makes it clear that whether a per- 
turbation lies in B or its complement depends not only on the amplitude 
of the perturbation but also - sensitively ~ on its direction in phase space. 
Moreover, even for perturbations in the right direction to leave B, a large 
amplitude is not necessarily more effective than a smaller one. 

I wish to thank Carlo Cossu for generously sharing his data with me. 
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